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The Bose-Einstein condensation of correlated atoms in a trap is studied by examining the effect 
of inter-particle correlations to one-body properties of atomic systems at zero temperature using a 
simplified formula for the correlated two body density distribution. Analytical expressions for the 
density distribution and rms radius of the atomic systems are derived using four different expres- 
sions of Jastrow type correlation function. In one case, in addition, the one-body density matrix, 
momentum distribution and kinetic energy are calculated analytically, while the natural orbitals 
and natural occupation numbers are also predicted in this case. Simple approximate expressions for 
the mean square radius and kinetic energy are also given. 
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I. INTRODUCTION 



The first theoretical prediction of the famous phenomenon known as Bose-Einstein condensation (BEC) was made 
in 1924 and 1925 by Bose ||l| and Einstein respectively. In a system of particles obeying Bose-Einstein statistics 
where the total number of particles is conserved, there should be a temperature below which a finite fraction of all the 
particles condense into the same one-particle state Seventy years later, in a remarkable experiment, Anderson et 

al. 1^ have cooled magnetically trapped ®^Rb gas to nanokelvin temperatures, and observed the BEC. This discovery 
has generated a huge amount of theoretical investigations ]l0|-|27t . 

The main feature of the trapped alkali-metal and atomic hydrogen systems (which obey the Bose-Einstein statistics) 
is that they are dilute, i.e., the average distance among the atoms is much larger than the range of the interaction. 
The crucial parameter defining the condition of diluteness is the gas parameter x = '^^'^ i where n is the density of the 
system and a is the s-wave scattering length [|l5| . There are two ways to bring x outside the regime of validity of the 
mean field description. The first one consists by increasing the density, while the second one consists by changing the 
effective size of the atoms [|l^ . In the present work we follow the second way. 

The characteristic dimension of a harmonic oscillator (HO) trap for s'^Rb is 5 = = (1 - 2) x lO"*! while 

the scattering length lies in the range 85ao < a < 140ao, where = 0.5292A is the Bohr radius. The atomic density 

in the trap is n ~ 10^^ — 10^** atoms/cm'^ giving an inter-atomic distance I = — lO'^A [Q. In this case, the 

effective atomic size is small compared both to the trap size and to the inter-atomic distance ensuring the diluteness 
of the gas. However, the effects of inter-particle interactions are of fundamental importance in the study of the BEC 
dilute-gas where the physics should be dominated by two-body collisions described in terms of the s-wave scattering 
length a. In the case of positive a, it is equivalent to consider a very dilute (atomic) system of hard spheres, whose 
diameter coincides with the scattering length itself . The natural starting point for studying the behavior of those 
systems is the theory of weakly interacting bosons which, for inhomogeneous systems, takes the form of the Gross- 
Pitaevskii theory p8| , p9| . This is a mean-field approach for the order parameter associated with the condensate. It 
provides closed and relatively simple equations for describing the relevant phenomena associated with BEC . 

In the present work, we avoid to use the Gross-Pitaevskii's theory and we study the BEC in a phenomenological 
way where the bose-gas is considered as a many body system [^^. In particular, we study the ground state of a 
system of correlated bosonic atoms at zero temperature, trapped by a HO potential. The key quantity for this effort 
is the two body density distribution (TBDD) p(ri,r2), which expresses the joint probability of finding two atoms 
at the positions ri and Y2, respectively. In the mean field case, the TBDD is the product p{vi,Y2) — p(ri)p(r2) 
where p(r) is the density distribution (DD) of the system. As the mean field approach fails to incorporate the inter- 
particle interactions which are necessary for the description of the correlated Bose system, we introduce the repulsive 
interactions through the Jastrow correlation function /(| ri — r2 |). In this approximation the TBDD has the form 
p(ri,r2) = N p{yi) p{y2) P {ri2) , where the normalization constant N ensures that / p(ri,r2)dridr2 = 1 [|l|,|2|j3|. 

The calculations of the ground state properties of the correlated Bose system (DD, rms radii, momentum distribution 
(MD) and kinetic energy (KE)) are made using four different expressions for the correlation function /(ri2). Two 
of them are Jastrow type Gaussian functions similar to those which have been used extensively in Nuclear Physics 
p2| , ^ , while, three of them are of hard sphere form and they are more realistic in describing atomic and molecular 
systems [^^p6|. The various expressions of the correlation function have one or two parameters (the correlation 
parameter (3 and the hard sphere radius rp). In the case of the Jastrow type Gaussian function, analytical expression 
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of the one-body density matrix (OBDM) is derived and is also used for the calculation of the corresponding natural 
orbitals (NO's) and natural occupation numbers (NON's). Finally, we give analytical or approximate expressions for 
the rms radius and the mean kinetic energy (T) (in one case) which depend on the parameters f3 or/and rg. 

The plan of the paper is as follows. In Sec. II the general definitions related to the Bose system are considered. 
Details of the lowest-order cluster expansion and analytical expressions of the DD (and in one case of the MD) are 
given in Sec. III. Numerical results are reported and discussed in Sec. IV, while the summary of the present work is 
given in Sec. V. 



II. GENERAL DEFINITIONS 



Let 'I'(ri, r2, • • • , r^) be the wave function which describes the inhomogeneous atomic system. In the case where 
this system is composed of bosonic atoms at zero temperature, all atoms occupy the same single-particle ground 
state. The many body ground state wave function is then a product of A identical single particle ground state wave 
functions. This ground state wave function is therefore called the condensate wave function or macroscopic wave 
function and has the form 

5'(ri,r2, ■■■ ^ta) = V'o(ri)'0o(r2) • • • V'o(rA), (1) 

where ipoij) is the normalized to 1 ground state single particle wave function describing a bosonic atom. It is worth to 
indicate that Eq. ([^) is valid even when weak interactions are included. In this case the wave function ^'(ri, r2, • • • , r^) 
is still, to a very good approximation, a product of A single particle wave functions which are now obtained from the 
solution of a non-linear Schrodinger equation, the well known Gross-Pitaevskii equation. 

The density distribution p(r), which can be directly observed in a non-destructive way, is defined as js^ 

p(r) = J «'*(r, r2, • • • , rA)*(r, r2, • • • , rA)dr2dr3 • • • dr^. (2) 

In the condensation case, where the wave function of the system is given by Eq. (|l]), the DD, using relation (^), takes 
the form 

p{r) =1 Mr) P . (3) 
The TBDD p(ri,r2), which is a key quantity in this work, is defined as 

p(ri,r2)= y**(ri,r2,---,rA)^'(ri,r2,---,rA)dr3---drA, (4) 

which in the condensation case, takes the simple form 

p(ri, r2) =1 Vo(ri) H Mr2) P= p(ri)p(r2). (5) 

The TBDD is needed for the evaluation of the two-body properties of the atomic system, mainly for the interaction 
energy. It gives also direct information about the correlations between the atoms of the system. This is the reason 
that the calculation of this quantity is very significant for the study of atomic or in general quantum many body 
systems. 

Another quantity characterizing the atomic system is the OBDM p{ri,r'i) l22U3a-Bl| which is defined as 



p(ri,ri) = y"**(ri,r2,-- •,rA)*(ri,r2,- • •,r^)dr2dr3- • - drA. (6) 

The knowledge of the OBDM is also very important because this quantity is connected to the position and to the 
momentum space. The diagonal part of the OBDM gives the DD p{r,r') = p{r), while the MD n(k) ||l0| , ^ , p3| is 
given by a particular Fourier transform of it 

"(k) = J p{ri,r[) exp [ik(ri - r'^)] dridr[. (7) 

In addition, the OBDM may be expanded in terms of its eigenfunctions ipiir) corresponding to the eigenvalues rii, 
that is 
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p(ri,r;) = ^n,^*(ri)V'.(r'i), (8) 

i 

where 

^n, = l. (9) 

i 

The eigenf unctions ipiir), which are called natural orbitals, and the eigenvalues rii, which are called natural occupation 
numbers, are obtained by diagonalizing the OBDM through the eigenvalue equation 

J p{rur[)Mr[)dr[ = n,?A,(ri). (10) 

The condition, generally adopted, for the existence of condensation is that there should be one eigenvalue which is 
of the order of the number of the particles in the trap. 



A. Harmonic oscillator trap and Gross-Pitaevskii equation 



The first step for the study of the BEC at zero temperature is to neglect the atom-atom interaction. In this case 
(named independent particle model (IPM)), we consider that the atoms are confined in an isotropic HO well and the 
Schrodinger equation takes the form 



2m 2 



V'(r) = £V'(r)> 



where the ground state single-particle wave function ?/'(r) has the form 



i'oir) 





■ r2" 
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(11) 



(12) 



and the width b is 



b = 



1/2 



The normalization of the wave function is J \ ip{r) dr — 1, while the DD has the form p(r) =| "00(1") p. 

It is obvious that for a system of non-interacting bosons in a HO trap the condensate has the Gaussian form of 
average width b. If the atoms are interacting, the shape of the condensate can be changed significantly with respect 
to the Gaussian 0|. The ground-state properties of the condensate, for weakly interacting atoms, are explained quite 
successfully by the non-linear equation, known as Gross-Pitaevskii equation, of the form 



2ni 2 m 



(13) 



where A is the number of the atoms, a is the scattering length of the interaction and /i is the chemical potential [Q. 
This equation has the form of a non-linear stationary Schrodinger equation, and it has been solved for several types 
of traps using different numerical methods |4^-|45[|. The presence of the third term, which is linear in A is responsible 
for the dependence of the gas parameter x on the density of the system. That dependence is negligible in the case of 
the IPM. In the phenomenological approach which we consider in this work, there is not a direct dependence between 
the condensation and the number of the atoms. The inter-particle correlations is incorporated in the mean field only 
by the correlation function which, in a way, depends on the effective size of the atoms. 



III. LOWEST-ORDER CLUSTER EXPANSION 



A dilute BE atomic system can be studied using the lowest-order approximation (LWOA) jl^]. In the LWOA the 
two-body density matrix (TBDM) has the form 
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p(ri,r2;r;,r2) = iVp(ri,r;)p(r2)/(|r;-r2 |)/(| ri - |), (14) 

where /(ri, is the Jastrow correlation function, which depends on the inter-particle distance and N is the normal- 
ization factor. 

The diagonal part of p(ri, r2; r'^, r2) is the TBDD 

p(ri,r2)=7Vp(ri)p(r2)/2(ri2), (15) 

while the OBDM is 

P(i"i,ri)=y" p(ri,r2;ri,r2)dr2. (16) 
The DD, which is the diagonal part of p{ri, r'l), can also be obtained from the integral 

p(r) f p(v)p{r2)f{\ r - r2 \)dr2 (17) 



In the mean field case, using Eq. (||), the OBDM becomes 

p(ri,r;)=Vo(ri)V'o(r'i), (18) 

where V'o(r) is the uncorrelated ground state wave function. In the case of the inclusion of the inter-particle interactions 
between the atoms, which give rise to the depletion of the condensate, the OBDM takes the form |Q 

p(ri,r'i) =noVS(ri)V'o(r'i) + ^n,^*(ri)^^(r'i), (19) 

where the sum X^i^^o "iV'i*(ri)'0i(ri) is the contribution arising from the atoms out of the condensate. 

The NO's V'i(ri) and the NON's ni are obtained by diagonalizing the OBDM through the eigenvalue Eq. (|l 
This is made expanding the OBDM in Legendre polynomials 

oo 

p{r,r') = p{r,r' , cos uj rr') pi{r,r')Pi{cosujrr'), (20) 

where 

21 + 1 

Pi(^i^') — 7, / p{r,r' , cos LOrr') Pl{cOSLUrr') d{cOSLUrr')i (21) 



and substituting Eq. ( pO| ) into Eq. (|T^) . Using the addition theorem of the spherical harmonics and integrating over 
the angle il^^ , the eigenvalue problem takes the form 

POO 

An / pi{r, r')y.Tir'yW = n^^^^^^rir), (22) 



where ifi^f ('') radial part of ipi (r) 

Mr) = ^Tir)Yim{^r) ■ 

In the present work, we use four different expressions of the correlation function f{r 12) for the study of the BE 
atomic systems. These expressions are 

Case 1 /(r) = 1 - exp [-f3r^] , (23) 
Case 2 /W = { " ' ; J :V (24) 
Cases /(.) = {'"'"P["o^^^'^'"''°^]' , (25) 



Case 4 f{r) = { ~o J ' ^ > ^0 ^ (26) 

0, r <ro 
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where r =| ri — r2 |. 

The correlation function /(r) in case 1 for large values of r goes to 1 and it goes to for r ^ 0. It is obvious that 
the effect of correlations introduced by the function /(r) becomes large when the correlation parameter (3 becomes 
small and vice versa. Case 2 is the classical hard sphere correlation function, which has been used extensively in 
atomic and molecular physics and depends on the hard sphere radius tq. Case 3, which is a generalization of case 1, 
depends also on the "hard sphere" radius rp of the interaction between the atoms. In the limit ro — + 0, /(r) goes to 
that of case 1. Finally, case 4 is a generalization of case 2 including a Gaussian function with the additional parameter 
(3. In the limit (3^0 f{r) goes to that of case 2. 

The above defined correlation functions were used to find analytical expressions of the DD through Eqs. (|l6|) or/and 
(|l7|). In case 1 the OBDM and the MD was calculated also analytically through Eq. (|l^) and (|^) respectivelly, while 
the NO'S and the NON's are calculated through Eq. (||). 



A. Analytical expressions in case 1 



In case 1 we found the analytical expressions of the DD and MD from the OBDM which have been calculated from 
Eq. ([l^ ) and has the form 



P(r,r') 



N 



7^3/2^3 

where N is the normalization factor of the form 



Oi(r, r') - 02i(r, r') - O^^iv, r') + 023(r, r')] 



N 



1 



1 



(l + 2y)3/2 (l + 4y)3/2_ 
and the one- and the two-body terms of the cluster expansion have the forms 



Oi(r,r') = exp 



02i(r,r') = 



(l+y)3/2 

022(r,r') = 02i(r',r), 



exp 



1 + 3j/ rl 
l + y 2 2 



(27) 
(28) 

(29) 

(30) 
(31) 



023(r,r') = 



(l + 2y)3/: 



■ exp 



-(1 + 2?/) 



exp 



1 + 2?/ 



(32) 



where ri, = r/b and y = j3b'^ . 

The analytical expression of the DD can be found from Eq. (p7|), putting r' — r, while the MD can be found 
analytically using Eq. (m). It takes the form 



n{k) 



exp 3^)3/2 exp 



1 + Sy 



1 



(l + 2y)3/2(l+4?;)3/- 



■ exp 



1 



l + 2y 



(33) 



where k}, — kb. 

The expressions of p{r) and n[k), given by Eq. ( p7| ) (for r' = r) and (|3^), respectively, have been used to find the 
analytical expressions of the mean square radius and kinetic energy of the trapped gas, through the integrals 



(T) = — 47r / n{k)k^dk, 
2m 



(34) 
(35) 



respectivelly. The expressions we found, for (r^) and (T), are 
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(T) = Nhuj 



-3 



1 + y 3 l + 2y 
(1 + 2j/)5/2 ^2(1 + 4y)5/2 



3 3 1 + 3y 

4 ~ 2 (l + 2y)5/2 ' 4(l + 4y)3/2 



3 l + 2y 

+ T- 



(36) 



(37) 



These expressions, which for a given HO trap arc functions of the correlation parameter y, could be used to define 
y from Eq. (|3^), if the rms radius of the trapped atoms is known and then to define the (T) from Eq (^) and vice 
versa. For very large values of y Eqs. (|3^) and ( ^ ) give the HO expressions of (r^) and (T), i.e. (r^) = and 
(T) — ^hu), respectively. 



B. Analytical expressions in case 2 



In case 2, where the correlation function is given by Eq. (|2J), the DD defined by Eq. (|17l), is written 

^ Oi(r) - ^ exp [-rg] ( Erf(rb + rofc) - Erf(-rb + r^b)) 



p{r) 



where N is the normalization factor of the form 



TV : 



1 - Erf 



V2 



and 



The one-body term Oi (r) has the form 

1 



Erf(z) = 



"* dt. 



Oi(r) = exp [-rl] 



1 - 2r^^ + ^ [~*^^'' ~ ^0^)^]) ^ 2 (^''^(''^ + ''0'') + Erf(-rb + rob)) 



(38) 



(39) 



(40) 



C. Analytical expressions in case 3 



In case 3, where the correlation function is given by Eq. (|2^), the DD defined by Eq. (|17|), is written 

N 



p{r) 



[Oi(r)-202(r,/3)+02(r,2/3)], 



537^3/2 

where N is the normalization factor of the form 

N ^[Ni~2N2{P) + N2{2p)r\ 

where 



Ni = l- Erf 



V2 



rob exp 



'o& 



exp [yrt^bl 
(1 + 2yfn 



1 - Erf rob 



l + 2y 



2 rob 



TT 1 + 2?; 



exp 



' Oh 



(41) 
(42) 

(43) 
(44) 



where rgb — rQjh. The factor N2(2(3) can be found from the factor N2{P) replacing y 2y. 

The one-body term is the same as in case 2 and is given by Eq. (Efl), while the two-body term 02{r,j3) has the 
form 



= (1 + ^)3/2 

exp[-rg] 



1 + 22/ 



2V V (1 + y)^/' 



1 + y 

[exp [-(rb -f rob)^] - exp [-{n - fofc)^ 



Erf 



-n + ?'ofc(i + y) 

(l + y)l/2 



2rb^/^(l -h 2/) 

The term 02{r, 2/3) can be found from the term 02{r, (3) replacing y —f 2y. 



(45) 
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D. Analytical expressions in case 4 



In case 4, where the correlation function is given by Eq. (|2^), the DD defined by Eq. (|17|), is written 

N 



p{r) 



^3^3/2 

where N is the normahzation factor of the form 



[Oi(r)-02(r,/3)], 



TV : 



1 - Erf 



(roh/V2) 



23/2 



yrob 



exp [-rlJ2] 



x/^ l + 2y 

The term Oi{r) is the same as in cases 2 and 3, while the two-body term is of the form 



Erf 



(l + y)i/2 



-Erf 



-Th + ?-ob(l + y) 
(l + y)V2 



(46) 



(47) 



(48) 



The above expressions of the DD (and of the MD in case 1) depend on two parameters in cases 1 and 2 and on three 
parameters in cases 2 and 4, i.e. the HO parameter &, the correlation parameter /? or/and the hard sphere radius tq. 
We may note, however that their dependence on /3 and is only through the dimensionless quantities y = /36^ and 
rob = ro/b. 



IV. RESULTS AND DISCUSSION 

The calculation of the DD of trapped Bose gas, confined in an isotropic HO potential with length b — lO'^A, has 



been carried out on the basis of Eqs. (|27|), (p8|), (41) and ( |46| ) which have been found using four different expressions 
for the correlation function. The dependence of the DD on the parameters y and tq are shown in Figs. 1 to 3. 

The DD, in case 1, for various values of the parameter y, including also the uncorrelated case (y = oo), has been 
plotted in Fig. la. It is seen that, the large values of y (y > 10) correspond to the Gaussian distribution (HO case), 
while when y becomes small enough (y < 1) the DD spreads out as in Gross-Pitaevskii's theory. For y > 10 the effect 
of correlations is small, while for very large correlations (y ^ 0.1) the DD is modified entirely from the Gaussian form 
which originates from the HO trap. The same effect of the inter-particle correlations on the DD appears also in case 
2 where the correlation function is that of the hard sphere. This is seen in Fig. lb where the DD has been plotted for 
various values of the hard sphere radius tq. The effect of the correlations becomes significant for tq > lOOOA. Thus 
the inter-particle correlations become large when tq approaches the characteristic trap length b. It is noted that, in 
Gross-Pitaevskii's theory, the difference between the correlated DD and the uncorrelated one becomes significant for 
smaller values of the parameter ro 0,11^ • 

In Fig. 2 the DD in case 3 for three values of the parameter y and various values of the parameter tq is shown. It is 
seen that small values of y result to a depression of the central part of the DD independent of the values of tq. This 
is obvious in Fig. 2, where for y = 0.1 and tq > 500A, the maximum of the DD is moved from the center [r/b = 0) 
to r/b ~ 1. The same behavior of the DD in case 4 can be seen in Fig. 3, where it is plotted for the same values of 
y and vq as in Fig. 2. The only difference is that the same depression of the DD in the central region appears for 
smaller values of the parameter y. 

The rms radius of the Bose gas has also been calculated analytically from Eq. (EbI) in case 1 and numerically in 



the other cases from Eq. (34) for various values of the parameters y or/and tq. The results are shown in Figs. 4 and 
5. From Fig. 4a, which corresponds to the DD calculated in case 1, it is seen that (r^)^^^ (vb — r/b) is a decreasing 
function of the parameter y and forw > 10 approaches the rms radius of the uncorrelated case. If we expand the 
expression of (r^)^/^, given by Eq. (pq), in powers of y~^ and truncate the expansion to the y~^ power, the following 
approximate expression is obtained, 

{rl) = 0.4366y"3/2 - 0.60432y-5/2 + 0.2541y'\ (49) 

The values of (r^)^/^ calculated from that expression are very close to the ones calculated from Eq. (^ for y > 2. 
The behavior of (r^)^^^, in case 2, as function of the parameter rp is shown in Fig 4b. In this case, (rf)^^^ is an 

increasing function of tq. We fitted the form (r^) = | -I- Cr^^^ to the values of the rms radius obtained in case 2 and 
we found the expression 
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{rl) = ^ + 0.3978r^/'. (50) 

This simple expression gives values for (r^)^/^ which are very close to the ones calculated from the DD as can be seen 
from Fig. 4b. 

The dependence of {rl)^^"^ on y and tq, in cases 3 and 4, are shown in Fig. 5, where it is plotted versus ro for three 
values of the parameter y. It is seen that, in both cases for the same value of y, (r^)^/^ is an increasing function of 
To, as in case 2, while for the same value of tq it is a decreasing function of y, as in case 1. This behavior of (r^)^/^ 
as well as Eqs. (^) and (|^ lead us to fit the numerical values of (r^)^/^ in cases 3 and 4 with two functions, which 
are combinations of the expressions of (r^)^/^ of cases 1 and 2. For the two cases we found the expressions 

Case 3 (rg) = ^ + {Q.2%l&y-^/'^ - 0.1290?/"^/2 ^. o.0323?/"^) + 0.1327^i]^^ (51) 

Case 4 (rg) = ^ + {hA^^Qy'^/'^ - 7.3299y-^/2 + 2.1458?/"^) rj^f . (52) 

From Fig. 5, it is seen that the numerical values of (r^)^/^ calculated in cases 3 and 4 from the DD are very close to 
the ones calculated from Eqs. ( ^l| ) and (|5|), respectively. 

The NO's and the NON's have been calculated, in case 1, by diagonalizing the OBDM through Eq. (^. The 
NON nis, gives directly the condensation fraction no as a result of the repulsive interaction between the atoms of the 
Bose gas at zero temperature. The NON's for the ground state and for some excited states are given in Table 1. It 
seems that, for strong correlations, a fraction of atoms spread out into many higher excited states. The condensation 
fraction no, versus the parameter l/y is plotted in Fig. 6. From that figure and from Table 1 it is seen that the 
effect of the correlations on no is small and all the atoms occupy the ground state, when y > 10. The effect of the 
correlations is prominent when y < 10, while the decrease of the parameter y (large correlations) induces a significant 
depletion of the condensated atoms exciting them into higher states. 

The NO's of the states Is, Ip and Id for y — 1 are shown in Fig. 7. It is seen that the interatomic correlations 
in the Is-state spread out the ground state wave function and consequently the condensation is appears in the outer 
region of the trap. From the same figure it is obvious that the NO's of excited states are much more localized in 
coordinate space than the equivalent HO orbitals. 

The MD in case 1 can be calculated analytically from Eq. ( |3^ ) or by Fourier transform of the NO's. The MD 
calculated analytically for various values of the parameter y has been plotted in Fig. 8a. It is seen that the large 
values oiy {y > 10) correspond to the Gaussian distribution, while when y becomes small enough {y < 1) the MD has 
a sharp maximum for fc = 0. The MD of the Is-state NO as well as of the residual excited states for y = 0.1 are shown 
and compared with the total MD in Fig. 8b. It is obvious that although the Is-state NO gives the main contribution 
to the MD, the NO's of the excited state contribute to the MD mainly in the large values of the momentum k. 

The dependence of the mean kinetic energy (T) on the parameters y calculated analytically in case 1, using Eq. 
(^), is presented in Fig. 9. It is seen that (T) has a maximum for y ~ 2.5. It is interesting to note that for the 
same value of the parameter y the NON's of the states Id and If have a maximum as can be seen from Table 1. The 
contribution of the Is-state NO and of all the excited states to (T) are shown in the same figure. It is seen that, for 
large values of the parameter y (weak correlations) the main contribution to (T) comes from the Is-state NO, while 
for large correlations there is a significant contribution coming from the NO's of the excited states. 

An approximate expression of (T), similar to the one of (r^) can be found if we expand the expression of (T), given 
by Eq. (|3^), in powers of and truncate the expansion to the y~^ power. The approximate expression has the 
form 

{T)/hLj = ^ + 0.1875?;^^/2 - 0.3355^"^/^ + 0.1092y"l (53) 

The values of (T) calculated from that expression are very close to the ones calculated from Eq. (^^ for y > 2.5. 

A final comment is appropriate. The Gross-Pitaevskii equation explains successfully the ground-state properties 
of the condensate for weakly interacting atoms, i.e. for small correlations between the atoms. The method of the 
present work can be used for larger correlations, provided that the expansions of the TBDM and the OBDM contain 
higher order terms. It should be noted also that, in the present work there is not a direct dependence between the 
condensation and the number of the atoms. The inter-particle correlations are incorporated in the mean field only 
by the correlation function which, in some way, depends on the effective size of the atoms. However, the condensate 
fraction can be found by calculating the natural orbitals and the natural occupation numbers by diagonalizing the 
OBDM from Eq. (|o|). These calculations have been made in case 1, while it is a very difficult task in the other cases. 
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V. SUMMARY 



The effect of ttie inter-particle correlations between Bose atoms at zero temperature is examined using a phenomeno- 
logical way to incorporate the atomic correlations. That is made by introducing the Jastrow correlation function in 
the TBDD. We examined the effect of correlations using four different correlation functions. The introduction of 
correlations change the shape of the DD and MD comparing with the Gaussian form, which corresponds to IPM. 
There is a reduction of the DD in the central region of the atomic system, and so an increase of the rms radius of the 
system, while the MD, which is calculated in case 1, increases in the region of small k and so there is a decrease of 
the mean kinetic energy of the system. In case 1 where the mean kinetic energy is calculated also analytically, there 
is a maximum for a certain value of the correlation parameter. Finally in case 1 the natural orbitals and the natural 
occupation numbers have been calculated and consequently the condensate fraction has been obtained for different 
values of the parameter y. 
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TABLE I. The values of the natural occupation numbers in case 1 for different values of the correlation parameter y . 



y nia nip nid ^i/ Sum 



100 0.99988 - - - 0.99988 

10 0.99634 0.00063 0.00042 0.00042 0.99781 

5 0.99055 0.00273 0.00108 0.00108 0.99544 

2.5 0.97771 0.00960 0.00186 0.00186 0.99103 

1 0.94422 0.03462 0.00172 0.00172 0.98228 

0.5 0.90815 0.06830 0.00082 0.00082 0.97809 

0.1 0.83097 0.15185 0.00001 0.00001 0.98284 

0.01 0.79273 0.19414 - - 0.98687 
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FIG. 1. The DD p(r) in cases 1 and 2 versus r for various values of the parameters y and ro, respectively, r is in units of 
the trap lenght 6, p{r) in and ro in A. The normalization is f p(r)dr = 1. 




FIG. 2. The DD in case 3 versus r for various values of the parameters y and ro. The units and the normalization are as in 
Fig. 1. 




FIG. 3. The DD in case 4 versus r for various values of the parameters y and ro. The units and the normalization are as in 
Fig. 1. 
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FIG. 4. The rms radius, (r^)^''^, in cases 1 and 2 versus the parameters y and ro, respectively. Case 1 corresponds to the 

analytical calculations of {r'^)^^^, while the points in case 2 correspond to the numerical calculations of it. The dashed line in 
case 2 corresponds to the fitting expression (50). {r^y^^ is in units of the trap lenght 6. 
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FIG. 5. The rms radii in cases 3 and 4 versus the parameter ro. The circles, squares and triangles correspond to the 
approximate calculations of (r^)^'''^ for three different values of the parameter y, while the dashed, dot and dashed-dot curves 
correspond to the fitting expressions (51) and (52) for the same values of y. (r^)^^^ is in units of the trap lenght h. 
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FIG. 6. The condensate fraction no, at zero temperature, for interacting atoms in case 1 versus the correlation parameter 




FIG. 7. The NO's (dashed hne) of the states Is, Ip and Id obtained by diagonahzation of the OBDM in case 1. The sohd 
hnes correspond to the HO wave-function with the trap length b. 
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FIG. 8. (a) The MD (normalized to 1) in case 1 versus k for various values of the parameters y. (b) The MD and the 
contribution to it from the NO's of the Is-state and of all the excited states for y = 0.1. k is in units of the inverse of the trap 
length (6"^) and the MD's are in A^. 
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FIG. 9. The mean kinetic energy per atom, (T), in case 1 versus the parameters y. The solid curve corresponds to the total 
values of (T), while the dashed line and the dotted line are the contributions to the total (T) of the NO's of the Is-state and 
of all the excited states, respectively. 
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